Sub-eV scalar dark matter through the super-renormalizable Higgs portal 
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The Higgs portal of the Standard Model provides the opportunity for coupling to a very light 
scalar field </> via the super-renormalizable operator (f>(H'' H). This allows for the existence of a very 
light scalar dark matter that has coherent interaction with the Standard Model particles and yet has 
its mass protected against radiative corrections. We analyze ensuing constraints from the fifth-force 
measurements, along with the cosmological requirements. We find that the detectable level of the 
fifth-force can be achieved in models with low infiationary scales, and certain amount of fine-tuning 
in the initial deviation of <j) from its minimum. 



I. INTRODUCTION 

About 95% of the energy budget of the Universe con- 
sists of "dark" - and unknown - components. This is 
a strong motivation for considering and studying hidden 
sectors beyond the Standard Model (SM). Gravitational 
effects of dark matter cannot reveal the mass of its con- 
stitutents, and indeed a wide variety of mass ranges, from 
the inverse galactic size to the super-Planckian scales, is 
conceivable. While many models that possess stable par- 
ticles with masses comparable to the SM energy scales 
have been a subject of incessant theoretical and experi- 
mental activity, models with light sub-eV mass scale dark 
matter received far less attention. 

Below the eV mass scale the dark matter would have 
to be of integer spin, and be produced non-thermally. 
The only chance of detecting such dark matter non- 
gravitationally would occur if such particles are converted 
into electromagnetic radiation in the external fields or 
they modify the interaction stength of SM particles. But 
if light dark matter interacts with the SM, then immedi- 
ately its lightness comes to question as the quantum loops 
with SM particle may easily destabilize the mass scale. A 
prominent particle in this category is the QCD axion [T] 
that interacts with the SM currents derivatively, j^d^a, 
and has its tiny mass generated by the non-perturbative 
QCD effects protected at any loop level. Because of the 
pseudoscalar nature of a and its derivative couplings, it 
does not generate a long-range attractive force. 

A very natural question to ask is whether SM allows 
for couplings to other types of sub-eV dark matter fields 
that lead to additional observable effects. For a recent 
review of the light sector phenomenology see, e.g. 
Real scalar field <f) and the vector field provide such 
opportunities with their couplings to the SM fields via 
the so-called Higgs and vector portals: 

( + \(f)'^)H^H Higgs portal ( 1 ) 



J^y^l.'^ dpi J II = Vector portal, 

where H is the Higgs doublet, A and A are parameters 
and is some locally conserved SM current, such as 
hypercharge of baryon current. If there is some initial 
value for (j) or fields with respect to their zero energy 
configurations, one can source part/all of the Universe's 
energy density from the coherent oscillations around the 
minimum. 

The perils of low mass scale stabilization are immedi- 
ately apparent in Eq. ([!]). Indeed, any loops of the SM 
fields would tend to induce the correction to the mass 
of (j) field ~ AA^y, where Aj/y is the highest energy 
scale in the problem serving as the ultra-violet cutoff. 
Therefore, A should be taken to incredibly small values, 
making this portal irrelevant for the phenomenology of 
sub-eV dark matter. In contrast, the vector portals and 
the super-renormalizable Higgs portal, AcjiH^ H, allow to 
avoid problems with technical naturallness. In the lat- 
ter case loop corrections scale only as log Kuv, while 
the quadratic divergences affect only the term linear in 
0, which can typically be absorbed in an overall field 
shift. In this paper we examine generic consequences of 
this coupling for the sub-eV scalar dark matter, leaving 
vector dark matter to future studies. 



II. SUPER-RENORMALIZABLE PORTAL TO 
THE SCALAR DARK MATTER 

The specific case of a singlet scalar (j) coupled via 
a super-renormalizable term of the type (j>H^H, (see 
e.g. [3HH] and references therein), has been mostly stud- 
ied in connection with electroweak and GeV-scale phe- 
nomenology, with a notable exception of [6l [9], where 
a possibility of super-weakly interacting Higgs-coupled 
dark matter was pointed out. The scalar potential in the 
model of interest reads as: 

V^--I^H^H + \iH^Hf+AH^H^+'^(t,\ (2) 



'Electronic address: fpiazza@perimctcrins titute.ca| 
t Electronic address: impospelov@perimeterinstitute.ca 



This model is explicitly renormalizable and does not re- 
quire any additional UV completion (if one is willing to 
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tolerate the usual fine-tuning problem with itself). 
We chose to redifine away possible linear terms in <j) by 
shifting the field, and absorbing A At/) into m\. 

After spontaneous symmetry breaking, the two fields 
acquire a vacuum expectation value, {H^H) = 
{(j)) = 00, where 



2\-A^/r 



2m^ 



(3) 



and V = 246 GeV. The potential ^ has a stable mini- 
mum only if /m^ < 2A, which is what we assume in the 
following; otherwise, it develops a runaway direction in 
the {4>, H) plane imless additional nonlinear 0** terms 
are introduced. The low energy dynamics is encoded in 
the two physical fields h and (p, defined as 



H 



1 







and with Lagrangian 



C = 



[dhf 



-{Av)hLp - 



A 



(4) 

(5) 
(6) 



As already noted, Higgs loops give only logarithmically 
divergent corrections to rriip. Therefore, the requirement 
of technical naturalness bounds the scale of m^p from be- 
low by the coupling A. In summary, by defining the 
dimensionless ratio x = A/m^, we assume x < 1 and 
X < a/2A, although also values x ^ I will be considered. 



In the above, gt/f are the Yukawa couplings to 
fermions. Those can either be fundamental, as the SM 
couplings to quarks and leptons, ghqq — rug/v, g^u = 
mi/v where niq {mi) is the mass of the quark (lepton) 
under consideration, or effective, as in the case of the 
nuclcons. The latter includes the contributions from all 
heavy quarks contributing to the coupling to gluons g^gg 
that provide a dominant contribution in the chiral limit 
[TU] . Below the QCD scale, the estimate of the effective 
Yukawa coupling of the Higgs to nucleons is rather un- 
certain due to a poorly known strangeness content of the 
nuclcon in the 0^ channel: 



200 



ghNN 



500 MeV 

V 



(9) 



This is much larger than the naive contribution of up and 
down quarks. 

The violation of EP is evident from the fact that the 
electrons and nucleons have couplings to the ip field that 
do not scale exactly with masses. 



g/te. 
me 



ghNN 



m„ 



(10) 



The effective coupling of the Higgs to the electromag- 
netic field, ghjj, is obtained by integrating out heavy 
charged particles, and the question of which one is 
"heavy" depends on the characteristic of (virtual) pho- 
tons. The coupling gh-yy can be written in the following 
form (see, e.g. ,12]): 




(11) 



III. FIFTH FORCE AND EQUIVALENCE 
PRINCIPLE VIOLATION 



The singlet (p couples to SM particles through the mix- 
ing with the Higgs field. Depending on the mass m^p and 
coupling A, the iy9-mediated attractive force can produce 
testable deviations from l/r^-gravitational force as well 
as composition dependence, thus violating the Equiva- 
lence Principle (EP). The leading contributions to p- 
couplings mediated by the cp-Higgs propagator is shown 
in Fig. [T] As a rule of thumb, the </3-couplings are sup- 
pressed with respect to the Higgs couplings by a factor 
of Av/mj^: 



Av 



gipx 



vghx 



(7) 



where guxx is the effective dimensionless coupling of 
the Higgs to x-particle at very low momentum tranfer. 
Therefore, the effective Lagrangian describing the inter- 
actions with the SM gauge and fermion fields takes the 
following form: 



off 



— [9hffJJ 
mt V 



(8) 



where summation goes over the quark and lepton fields 
with charges Qq and Qi, and the last term is due to 
the the M^-bosons. For the purpose of calculating the 
77 decay, one has to sum over e, ji, r and c, 6, t. 
Corrections coming from the light quark sector are sub- 
dominant, because in the chiral limit they contribute at 
two-loops. In practice, their contribution would amount 



h 



-->< Qh 



Av m 



h 




Figure 1: The mixing with the Higgs Av mediates the cou- 
pling of (f to SM particles. 
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to at most 10% correction. Including these fermion con- 
tributions gives gh-yyil^ = w^) ~ Q;EM/(87r). For the 
purpose of calculating the coupling of (p to nuclei when 
the EM fraction of energy is taken into account, electrons 
should not be included in the sum, and muon contribu- 
tion should include a form-factor. We are not going to 
pursue this calculation, because it turns out that gh-f-y 
provides a subleading contribution to the EP violation. 

Field if mediates a fifth force of range ~ iti~^. More 
precisely, at the Newtonian level of approximation, the 
total effective gravitational potential between two bodies 
A and B at relative distance r, presents a Yukawa con- 
tribution due to the interaction of the long range field 



V{r) 



niAmB , r\ 



(12) 



The scalar couplings a can be expressed in terms of the 
log-derivative of the masses as 



OlA 



V2Mp 



din iriAi^) 
dip ' 



(13) 



where Mp is the reduced Planck mass and mAip') in- 
cludes terms in the Lagrangian that are bilinear in the 
fields and couple to ip, such as those in eq. ([s]). When 
calculating a a, one should consider the leading univer- 
sal contribution from the nucleons and all the corrections 
that are specific to the element A (See e.g. [13] )• The 
main, species-independent part of the nuclear mass is 
given by ?7inuc(-^A + Za), and the universal coupling a is 
obtained from eqs. (|8|, ^ and ([13]): 



ghNN - 



IQ- 



V2Mp Av 

2 

^nuc '^h 

mh 
115 GeV 



(14) 



A 



10-8cV 



In the limit of a very long range force, the value of 
a is bounded by post-Newtonian tests of General Rela- 
tivity to < 10~^ [Hj. However, one can easily see 
that for mass range of below 10~^^ eV, the rela- 
tive strength of the ^-induced force drops below 10"^'' 
from the gravitational field strength, which would make 
it extremely challenging for experimental detection and 
immune to the Solar System tests. Thus, it is more in- 
teresting to consider intermediate-range forces. Tests of 
gravitational inverse-square law limit the Yukawa com- 
ponent of the gravitational potential [HI [16] . By means 
of equation ( [T4| , such tests give a bound on A. This 
is shown in Fig. [2] The two panels are elaborations of 
plots taken from Refs. [T3] and [T^]. A force with similar 
values of and ^ (a; ~ 1) is excluded in the range of 
masses ~ lO^^eV — lO^'^ eV. 

The calculations of the EP-violating part of the scalar 
exchange is a far more delicate excercise. One should rec- 
ognize that the equivalence principle is violated already 
at the level of nucleons, that is (7;i„„/m„ ^ ghpp/Tnp- As 
is well-known, the neutron and proton mass difference 



1 — V — I — I — I — I — 
excluded region 




Figure 2: We plot the constraints on the mass and cou- 
pling A — xm^p coming from fifth force experiments, and tak- 
ing ghNN to the maximum of its allowed range. The range 
of the force is just A = m^^ . The coupling a is obtained 
in eq. ( 14 1 by assuming ~ 120 GeV. For two diflterent 
mass ranges, the lines corresponding to x — 1, x = 10~^ and 
X — 10"'' are superimposed on the plots of references [15j 
(upper panel) and [TB] (lower panel). 



comes about because of the unequal quark masses, and 
electromagnetic contribution to the nucleon mass. One 
can estimate (m„ — 'mp)\m^^mci — 2.1 MeV and (r7i„ — 
TOp)|em ~ —0.8 MeV, so that together both contributions 
combine to the observable mass difference Am„p = 1.3 
MeV. The (ys-dependence of both pieces is completely dif- 
ferent. Because of the loop smallness of gh^yy the electro- 
magnetic fraction of nucleon mass is far less dependent on 
ip: d{mn-mp)\EM/dh ^ d{mn-mp)\m^^ma/dh. There- 
fore, when we estimate the mass of an atom, we add to 
the universal term proportional to the baryon number a 
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correction proportional the the nucleon mass difference: the field from its minimum at the time when ~ 3H. 



N - Z 

{N + Z)mnuc{^) H z — Anpni{ip) 



(15) 



The first term in (15) produces the universal coupling a 
calculated in (14 1. The composition-dependent correc- 
tion reads 



EPV ^ 



N — Z Am„j; 
2{N + Z) rriN 



mjv dAninp/dh 



9hNN 



Arrir, 



2{N + Z) 



X 3 X 10^3(16) 



This may lead to a sizable variation of acceleration Aa 
between light atoms with Z = N and heavy atoms with 

N-Z _ Q . 
2{N+Z) — 



Aa 



(17) 



Other important effect should be related to the depen- 
dence of the nuclear binding energy on ip, that can easily 
reach a level comparable to (16). More detailed consid- 



erations of nuclear mass dependence on Lp go outside the 
scope of the present paper. 

As long as we adhere to our naturalness condition 
A ~ m^, the present bounds on composition depen- 
dent EP violations {Aa/a < 10"^"^) are easily evaded. 
When the Earth is the common attractor of the two free- 



falling bodies, the relevant range m< 



10 km turns 



into extremely tiny values for the coupling A. Still, if we 
were to consider more fine-tuned scenarios (m^ <^ A), 
it is interesting to note that a fifth force attached to the 
Higgs portal displays a peculiar relation between com- 
position independent and composition dependent effects, 
as clearly follows from eq. (16). In principle, this al- 
lows to distinguish between the Higgs portal and, e.g., 
the string- inspired scenarios [131 UHl HH] • 



IV. COSMOLOGICAL CONSTRAINTS 

Since all couplings to SM particles are suppressed by 
a factor Av/mf^, the scalar field is sufficiently stable to 
be a non-thermal relic (for earlier studies of scalar dark 
matter see, e.g., [inilll])- Its decay rate into photons is 
smaller than the current Hubble rate as long as the mass 
is under a keV: 



,3 42 2 

47rm^ 



10--eVx.^(^)7M^V 

(18) 

However, as emphasized in [5], the constraints from the 
gamma ray background would provide much tighter con- 
straints, and in what follows we will concentrate on the 
sub-eV range. 

The abundance of (/j-particles today can be esti- 
mated [221 123| in terms of the initial misalignment ip^ of 



At that moment, the field starts oscillating around the 
minimum of its potential and behaves like non-relativistic 
matter. The number of particles in a comoving volume is 
conserved, so that n^/s — const, where s — O.A4:g^,T^ is 
the entropy density, is the number density of ip parti- 
cles and g.^, is the number of effective degrees of freedom in 
equilibrium with the photons. The (average) energy den- 
sity of p today is thus given by p° = m^pn^s'^ /s, where 
and s should be taken at t^. Unsing the relations 



.^..^ = mlpl/2 and s°/s, = {2/ g,){T° /T^f together 
with the Hubble rate 



Tj -1/2 

H ~ g,' 



T2 



(19) 



3Mp 

we express in terms of the parameters of the model to 
obtain 



-3/4 

9* 



1/2 



where is the number of degrees of freedom relevant 
for estimating H. Since T, ~ 10^{m^/eVy^'^GeV, in the 

mass range of interest — and O.Ag^^"^ /g^ ~ 0(0. 1). 

An important constraint on the model comes from the 
smallest allowed mass for a dark matter particle. The 
observations of smallest halos show their size to be com- 
parable to 1 kpc [24], which means that the Compton 
wavelength of p field would have to be comparable or 
smaller that this scale. This in turn imposes the con- 
straint on pi^,: 

1/2 

(21) 



<2x IOI8 GeVx 



0.1 



Notice that at the boundary of the allowed value, m^p ~ 
10~^^ eV, the oscillations start around 10 eV, that is just 
before the matter-radiation equality. 

By eq. ( [20| it is clear that the abundance of (/s-particles 
depends on the VEV of the field at the moment when the 
Hubble parameter becomes of the order of its mass. This 
is ultimately a matter of initial conditions. However, 
it is interesting to study the preceding evolution of p 
up to electroweak symmetry breaking. During inflation, 
while the field would classically stay constant, its vacuum 
expectation value gets random kicks of order H/2'k every 
Hubble time due to quantum fluctuations (see e.g. [25J). 
Its behavior can be described formally with a Langevin 
type equation [TB] 



dp 



Hip) 
2tt 



C(P). 



(22) 



In the above p = In a is the number of c-folds and ^ is a 
Gaussian random variable. Its p-averages are {C{p)) — 0, 
{C{p)C{p')) — S{p — p'). It is straightforward to estimate 
the expected shift in the field during inflation: 



|A^i„f| = v/((0c„d - Cj^inr) = ^ (^J^"^'dp'H^p') 



1/2 



(23) 
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As a rough order of magnitude, this gives Ac/ijuf > 
10i?cMB, where Hqmb is the Hubble parameter at the 
epoch where the scales relevant for the CMB left the 
horizon. However, in scenarios with a long epoch of self- 
regenerating inflation, Acjjinf can be much larger. 

At the onset of radiation domination quantum fluctu- 
ations become irrelevant and the field is governed by the 
classical equation 



3H4) 



A{H^H) 0. 



(24) 



The behavior of cj) up to electro-weak phase transition 
is obtained by neglecting the second derivative and the 
mass term from the above equation. While the Higgs 
field is in thermal equilibrium we have, with good ap- 
proximation [23], (WH) = 3T^. By using (191 we thus 
get 



'3AMp . 



(25) 



To a good approximation, the field has a constant veloc- 
ity, which justifies neglecting the second derivative term 
in ( 24 ) . Every time a relativistic species leaves the ther- 



mal bath, 5* decreases, giving a little "kick" to the field's 
velocity. Thus, the details of this mechanism depend on 
the physics beyond the Standard Model. By making the 
minimal assumption ~ 100 we count only for the SM 
degrees of freedom and obtain a shift in field space. 



'-'EW — Vend 
Mp 



-0.4 AtEw ^-2 XX 10 



-6 "V 



lO-iOeV 



(26) 

where the subscript EW indicates quantities at electro- 
weak phase transition. 

Finally, at the onset of EW phase transition, the field 
finds itself displaced from its true minimum by an amount 
ipEw = 4>EW — 00, where 0o is given in eq. ([3|, 

^°"-^^^^'^"(io^)"^^^- 



A potential disparity between (21 ) and (27) signifies pos- 



sible fine-tuning problem. The starting point for the 
field at the end of inflation would have to be reasonably 
close to 4>Q, Therefore, if we start, say, with t^end = at 
the end of inflation, the field starts running towards its 
true minimum, 00 thanks to the coupling to the Higgs 
(26). However, for masses < 10~^ eV, the shift dur- 



ing radiation domination ( 26 ) is irrelevant with respect 



to the scale set by (27). If Lp mediates long/intermediate 
range (A > cm) forces, its initial value has to be fine 
tuned ((,5* <C |0o|), or otherwise (yS-particles are overpro- 
duced. 

Any super-cold dark matter, such as axion or field 
discussed in this paper, are prone to the CMB constraints 
on the amount of isocurvature perturbations (For the re- 
cent discussions of the axion isocurvature perturbations 
see e.g. [HI [57]). Before going to implications of these 
constraints for the model, we would like to comment that 



the scalar field with the quadratic potential is less sus- 
ceptible to the isocurvature constraints than axions. In 
case of the quadratic potential the increase in the homo- 
geneous displacement from the minimum, ^p.^,, over the 
fluctuating value 5<j) leads to the Slp/lp^, suppression of the 
isocurvature perturbations, that in principle can be made 
arbitrarily small by the increase of In contrast, the 
increase in the homogenous value of the axion field due 
to the periodicity of the potential Va[a) = V{a + 27^ fa) 
can lead to at most 5a/ fa suppression of the isocurvature 
perturbations. 

During infiation, the field undergoes fluctuations of 
order Sp — H/2'it as any other light held, H being the 
Hubble rate at the time when the fluctuation exits the 
horizon. The produced perturbation are of isocurvature 
type. Following the standard treatment that also applies 
to axions [531 12H] , we can estimate the power spectrum 
of entropy perturbations Vs{k) and compare it to that 
of curvature perturbations 'P-R,{k). The ratio of the two 
defines a parameter 



a{k) _ Vs{k) 



a{k) Vnik) 



5ec, 



(28) 



where f2c is proportional to the total energy density in 
dark matter, e is the usual inflationary slow-roll parame- 
ter and Cs the speed of sound of the adiabatic fluctuations 
during inflation [29]. 

By using (20) and ^Ich^ ~ 0.1 we get, for small a(A:o), 



.(.)=4.7xl0".c.2^(j53^) 



,1/2 



(29) 



The above result can nicely be re-expressed in terms 
of the tensor-to-scalar ratio r = 16eCs. At the 
pivot wavenumber ko — 0.002 Mpc~^ the limit set by 
WMAP-fBAO-fSN is a(fco) < 0.067 '28'. This gives the 
rather strict constraint 



fir 



( "V A 
VlO-9 eV/ 



1/2 



< 2.3 X 10 



-10 



(30) 



which is very similar to the conclusions reached for the 
axion cosmology [IHl 127] ■ If we insist on fl^p making most 
of the cold dark matter density, this result shows that 
the detectable level of inflationary gravitational waves 
(r > 10^^) implies a very light scalar, close to the bound 



(21), which would make mediated flfth force totally 
negligible. Conversely, a detectable level of the fifth force 
{m^ > A > 10^^ eV) would imply tiny r on the order of 
2Q-10 favoring some intermediate scale infiationary sce- 
narios, H ^ 0{r^^^ X 10^^) GeV. Given that in some 
models of inflation (see e.g. [30]) the Hubble parameter 
can be as low as H ^ GeV, producing a tensor to scalar 
ratio r ~ 10~^*, constructing an inflationary model with 
the Hubble parameter at some intermediate scale does 
not pause any model-building challenge. 
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V. DISCUSSION 

The model we considered in this work is very simi- 
lar to the linearized version of the Brans-Dicke (BD) 
theory when the scalar field is supplied with the mass 
term. Indeed, the transformation from the Jordan to 
the Einstein frame puts the BD scalar in front of any 
dimensionful parameter. Therefore, the A parameter 
from the model considered here can be identified with 
A ^ TO^/(a;^/^Mp), where uj is the BD parameter. It 
is very important to keep in mind, however, one cru- 
cial difference. In the BD theory, the (/)-field also cou- 
ples to all massive states that may exist beyond the SM 
states, and therefore, even at the electroweak scale one 
should expect the extension of Eq. ^ by additional 
higher-dimensional operators. Such terms alter the cou- 
plings of BD scalar to matter, and make couplings to 
gauge bosons, e.g. g^^^, different from the values in the 
model considered here. Moreover, the BD theory requires 
explicit UV completion, while the model with coupling 
via the super-renormalizable portal assumes that higher- 
dimensional operators are absent from the beginning and 
generated only via the SM loops with the (/^-independent 
UV cutoff. 

The key feature of the model considered here is its 
technical naturalness. It allows to have a relatively light 
scalar dark matter that generate medium-range attrac- 



tive force without extra fine tuning of the parameters 
in the Lagrangian. A detectable level of the fifth force, 
would have to be combined with inflationary scenarios 
with low r and face with the potential fine-tuning prob- 
lem in initial condition of the scalar field value. One 
of the most interesting (albeit fine-tuned) scenarios that 
can have particle physics implications not considered in 
this paper is the (/)-dependence of the electroweak phase 
transition. If is taken comparable to the Hubble rate 
at T = 100 GeV, A > can lead to - 10^ 

GeV^, thus altering the properties of the electorweak sec- 
tor close to the phase tansition point. This way, one could 
change the order of the phase transition, and make it first 
order if the effective Higgs mass is pushed below 50 GeV. 

The model considered here falls into the class of the 
"super-cool" dark matter models, such as axion dark 
matter. Another example, worth of investigation is the 
vector dark matter. There, the coupling of vector fields 
to the SM and the mass of the vector fields do not have to 
follow the strength xrange=const constraint of the scalar 
case. This could open more room for the fifth-force me- 
diated by the vector-like sub-eV dark matter. 

We would like to thank N. Afshordi, N. Barnaby, J. 
Bond, C. Burgess, A. Erickcek and A. Nicolis for use- 
ful discussions. Research at the Perimeter Institute sup- 
ported in part by the Government of Canada through 
NSERC and by the Province of Ontario through MEDT. 
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